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Abstract 



We show that the highly excited rovibrational spectra of a diatomic 
molecule and the closely related slow atomic collision processes contain more 
systematics and require less parameters to characterize than the Rydberg 
spectrum of an atom. In the case of a single channel, for example, we show 
that a single short-range parameter gives a complete description of slow col- 
lisions for practically all angular momenta, and covers an energy range of 
hundreds of millikelvins. The same parameter also describes the highly ex- 
cited rovibrational spectra in the threshold region, including states of different 
angular momenta. Sample applications and predictions of the theory are dis- 
cussed, including comparisons with experiment. 
PACS number (s): 34.10.+x,33.20.Vq,32.80.Pj 
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This paper examines the properties of a two-atom system in a small energy region around 
a dissociation threshold, with the region above corresponding to slow atomic collisions, and 
the region below containing the highly excited rovibrational states of the corresponding 
diatomic molecule. We discuss a type of systematics that relate states of different relative 
angular momenta, thus greatly simplifying our understanding of such systems. Comparison 
of the theory with existing experimental data [p]] also provides the first verification of the 
breakdown of the large-quantum-number formulation of the correspondence principle |2]-[1 • 

With recent developments such as the analytic solutions of the Schrodinger equation 
for 1/r 6 and 1/r 3 potentials PJ^J^l j our understanding of two-atom systems is already con- 
ceptually comparable to the quantum defect theory (QDT) of atomic Rydberg spectra 0. 
Namely, the slow atomic scattering and the rovibrational spectrum of a particular angular 
momentum can be understood in terms of the long range solutions and a parameter that is 
a slowly varying function of energy in the threshold region |5|,§,|)|,||f|] . In this work, we show 
that a two-atom system possesses an even greater degree of systematics in the following 
sense. Unlike the QDT for atomic Rydberg spectra in which different angular momentum 
states have different quantum defects with no general relationships among them, a single 
parameter is sufficient to characterize slow atomic collisions for practically all angular mo- 
menta. The same parameter also characterizes the rovibrational spectra in the threshold 
region, including states of different angular momenta. In other words, in addition to the 
relationship between the bound spectrum and scattering, as expected from traditional QDT 
formulations, a two-atom system has also the unique property that scattering of different 
angular momenta are related, and so are the bound spectra of different angular momenta. 
While a similar idea can be implemented in a numerical fashion , it becomes a much more 
powerful tool in the present formulation. 

The origin of this systematics is not difficult to understand and is due to a combination 
of the following three properties of a typical molecular system, (a) Atoms are strongly 
repulsive at short distances, (b) Atoms are heavy compared to electrons, (c) The atom- 
atom interaction at large distances behaves as 1/r 11 with n > 2. 
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The combination of the first two properties gives rise to the well-known separation of 
tightly bound rovibrational states into the product of two parts ip v ip r , in which the vibra- 
tional wave function ip v is, to the lowest order, independent of rotational quantum numbers. 
For rovibrational states that are highly excited, the same properties ensure that the radial 
wave function, up to a normalization constant, remain nearly independent of the angular 
momentum until a distance is reached where the rotational energy term becomes compara- 
ble to the electronic and other energy terms. This characteristic is uniquely molecular and 
does not apply to the electronic wave functions of an atom, for which the rotational term 
dominates the behavior at sufficiently small r. 

The implication of the third property can be understood from qualitative behaviors of 
the effective potential —C n /r n + h 2 l(l + l)/(2/xr 2 ), in which the first term represents the 
long range interaction and the second term represents the centrifugal barrier. For n > 2 and 
I 7^ 0, this potential has a maximum and crosses zero at r x = [1(1 + i)] _1 /( n-2 ) [3 n} where (3 n 
is a length scale associated with the —C n /r n interaction as defined by 

(3 n ee (2^C n /h 2 ) 1/(n - 2) ■ (1) 

From a physical point of view, r x specifies the location at which the rotational potential is 
equal to the long range interaction. For an / that is not too large, r x is of the order of f3 n . 
Since n is greater than two, for distances that are smaller than r x , the long range interaction 
quickly dominates, and correspondingly, the importance of the angular momentum quickly 
diminishes. Mathematically speaking, this means that for a potential that goes to zero 
faster than 1/r 2 at large distances, a pair of linearly independent solutions exist, which, 
at relatively small distances, are not only independent of energy, but also independent of 
angular momentum. 

The combination of these characteristics leads to the following important conclusion. 
With proper choice of long range solutions, a quantum defect theory for molecular rovibra- 
tional states and slow atomic collisions can be formulated in which the short-range param- 
eters are not only nearly independent of energy, but also nearly independent of the relative 
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angular momentum. This conclusion is very general and powerful. Its formulation and ap- 
plication are illustrated here for a single channel problem with a van der Waals asymptotic 
interaction. The formulation for attractive 1/r 3 interaction and multichannel generalizations 
will be presented elsewhere flIQ| . 

For a single channel with an asymptotic van der Waals interaction (— Cg/r 6 ), the proper 
pair of long range solutions are ones with the behavior 

f d ^ (2/vr) 1 / V/^r 1 / 2 cos Q(r/&)- 2 - |) , (2) 
& ^ -{2/,fl\r/^ sin Q(r//? 6 )- 2 - , (3) 

for both positive and negative energies. This pair, which has not only energy-independent, 
but also angular-momentum-independent behavior at small distances, is related to the f° 
and g° pair defined in || by a linear transformation. The fact that such a pair exist is a 
result of the exponent n that characterizes the long range interaction, —C n /r n , being greater 
than 2. 

With this choice of long range solutions, the single channel quantum defect theory is still 
formally the same as discussed previously |||6],[§]. 111 particular, the bound spectrum of any 
potential with the behavior of V(r) — > —Cq/t 6 at large distances is still given rigorously by 
the crossing points between a x function and a short-range K matrix ||{|] : 

(e.) = K c (e, I) . (4) 



Here e s is a scaled bound state energy defined by 

1 e 



(5) 



s 16(7i 2 /2//)(l//3 6 ) 2 ' 

The x c function, corresponding to the choice of long range solutions as specified by Eqs. (0) 
and (0), is given by 

c( v (Y d /X d ) + tan(7rz//2)(l + M d )/(1 - M d ) 
Xl[€s> l-(Y d /X d )tan(nu/2)(l + M d )/(l-M d ) ' W 
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in which M e i = G e i(—v)/G e i(u), with v and G e i being defined in 0. Plots of Xi functions 
for different I are shown in Figure |I[ K c is the short range K matrix that results from the 
matching of the short range solution and the long range solution. It is given specifically by 

KH n (f!l\ (f!l'/f!l) - (Ud/Ud) m 

{ , ) WJ (9d'/9 c d)-(u d '/u d ) ' [{) 

evaluated at any radius beyond r at which the potential has become well represented by 
-C n /r n . 

Above the threshold, the scattering K matrix is given by a similar expression as derived 
previously |J|,|: 

K x = tan 5, = (K c Z c gg - Z%)(ZJ f - K^Z^y 1 , (8) 

with the Z c matrix being related to the Z matrix defined in || by a linear transformation. 

The key difference between the present formulation and the previous one is the following. 
With the choice of long range solution pair as specified by Eqs. (0) and (||), and the property 
of the wave function discussed earlier, the parameter K c is, under the condition of (3 n ^> r 
not only independent of energy, but also independent of angular momentum / in 
the threshold region [see Eq. (|^)]. In other words, a single parameter, K c , can describe 
scattering and bound states of different angular momenta. 

The determination of the parameter K c is, in principle, straightforward theoretically. It 
can be done, for example, by solving the radial equation at a small energy (include zero 
energy) and matching it to the proper long range solution. However, since our present 
knowledge about both the short range interaction and the C n coefficients is still not suffi- 
ciently accurate for many systems, we focus here on the direct experimental determination 
of both the K c parameter and the C n coefficient. In particular, we show that if the C n 
coefficient that characterizes the dominant long range interaction is known accurately, the 
parameter K c can be obtained from the measurement of a single binding energy (similar to 
the determination of Kf discussed previously [§]). If binding energies of more than one state 
are known, in addition to the prediction of the parameter K c , an accurate determination of 
C n can also be made, especially when the two states are closely spaced in energy. 



In Table |L the first column presents the experimental results [l| for states of Rb2 



characterized by quantum numbers F\ = 3, F 2 = 3, F = 6, I = 2, and T = 8 flTIJI^]. The 
results in the column labeled by Theory II represents a simple calculation making use of 
only a single experimental energy. In this calculation, the energy of the level labeled by 
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1 in Table | is first scaled according to Eqs. ([5]) and ([I]) with /i = 77392.368 a.u. 



and a value of C@ = 4426 a.u. from the theoretical calculation of Marinescu et al. 0]. The 
parameter K c is determined by evaluating the Xi=2 function at this scaled energy ||, which 
yields a value of K c = 0.2841. The crossing points of this constant with xf =2 (e s ) give the 
predictions of the other energy levels listed in this column (see Fig. p]). Note that this simple 
calculation already yields excellent agreement with experiment and compares very favorably 
with a much more complex numerical calculation [|IJ listed under Theory I in Table |. This 
excellent agreement shows that the parameter K c is indeed, to a very good approximation, 
a constant in the energy range covered by the experiment (~ 10 GHz). This is further 
confirmed by the fact that enforcing K c to be constant leads to an even better agreement, 
a procedure that demonstrates a new method for the experimental determination of the 
leading C n coefficient from a minimum of only two binding energy measurements. 

Keeping in mind that the Cq value used in the first calculation is not necessarily the 
true C*6, it is allowed to vary in the following calculation. With more than one experimental 
binding energy available, we vary Cq in such a way that the two x c 's evaluated at the two 
different scaled energies (they may also correspond to different I.) are the same, in other 
words, we force the parameter K c to be a constant. The refined Cq is then given by the root 
of the following equation 

Xf(e.i)-X?'(<«) = 0, (9) 

where e s i and e S 2 are two experimental energies scaled according to Eqs. (||) and (0) for van 
der Waals interactions. The reduced masses are known with great precision from atomic 
masses so that the only unknown on the left hand side is the Cq coefficient. This procedure, 
when applied to the two d states labeled in Table | by v max — v = 1 and 2, respectively, leads to 



a revised C$ = 4533 a.u., which is in good agreement with Cq = 4550 ± 100 a.u., determined 



by Boesten et al. |nj . This C§ value refines the energy scaling and leads to a revised value 
of K c = 0.3106. The column in Table | labeled as Theory III gives the result of energy 
levels predicted using this set of revised parameters, and an even better agreement with 
experiment is achieved. This procedure demonstrates a new method for the determination 
of the C n coefficient from two or more values of experimental binding energies. It requires 
no knowledge of the short range interactions, nor does it require numerical solutions of the 
Schrodinger equation. The closer in energy those two states are, the better this method 
works. 

As stated earlier, the parameter K c describes not only the rovibrational states of a 
particular relative angular momentum, but also the other angular momentum states in the 
threshold region (see Fig. p. Table || gives the predictions of the bound spectra of 85 Rb2 
characterized by quantum numbers F\ = 3, F 2 = 3, F = 6,l,T = F + I, with I = 0, 2, 4, 6, 
respectively |TT|. They are predicted using parameters C 6 = 4533 a.u. and K c = 0.3106. 

The fact that K c is nearly a constant in the threshold region has also the following 
important implications, (a) It provides an experimental verification of the breakdown of 
the large-quantum-number formulation of the correspondence principle, as a semiclassical 
theory would have predicted different results for, e.g., the spacing between scaled energies 
||. (b) Since 1 GHz corresponds to 47.99 mK, the fact that K c is independent of energy 
over a range of 10 GHz below the threshold implies that the same constant K c can describe 
collisions over hundreds of millikelvins above the threshold (for all Vs that are not too 
large). (Similar conclusion can also be derived for Li |j,[nj). All the energy and angular 
momentum dependences in this energy range, including shape resonances, are determined by 
the dominant long range interaction, in combination with the centrifugal barrier. And they 
are described analytically by the solutions for the long range potentials [§]||(1|]. This result, 
along with the fact that for small energies, the scattering for large angular momentum 
is nearly independent of the value of K c and is determined entirely by the long range 
interaction, means also that the same parameter describes slow collisions for practically all 



angular momenta. As one implication, a single short-range parameter (a few parameters in 
the case of multichannels pUfl ) describes not only the ultracold region of energies, but also the 
region of evaporative cooling and the region of laser cooling. In contrast, in the traditional 
effective range expansion, a single parameter, the scattering length, gives little more than 
the cross section at zero energy. Figure |2| shows the partial scattering cross sections for s, d, 
and g waves predicted by the same short range parameter K c = 0.3106 for two spin-polarized 
85 Rb atoms in state F\ = 3, Mp% = 3, F2 = 3, and Mp2 = 3. Only a small energy range is 
shown here to make for easy identification of the narrow g wave shape resonance, observed 



previously by Boesten et al. |14j . More data on the scattering cross sections, including those 



for larger angular momenta and over a greater energy range, and the summation formulae for 
summing over all angular momentum states will be presented elsewhere [L0[|. (c) Since the 
hyperfine splitting is typically of the order of 1 GHz, it means that the frame-transformation 
method for atom-atom scattering [0, when properly formulated with K c (or the Kf 



being the short range parameter flOfl , should work well in a multichannel formulation that 
includes hyperfine structures [[Tall- 



in conclusion, a theory of slow atomic collisions and molecular rovibrational states has 
been formulated that takes full advantage of the intrinsic characteristics of a typical molec- 
ular system. This formulation establishes a simple relationship between states of different 
relative angular momenta and greatly reduces the number of parameters required to describe 
slow atomic collisions and highly excited molecular rovibrational spectra. Multichannel gen- 
eralizations do not change the key concepts in any substantial way, and the theory can also 
be implemented in exact numerical calculations. As far as parameterization is concerned, 
the theory can be further extended to an even greater range of energies by, e.g., taking into 
account the energy and angular momentum dependences induced by the correction terms 
to the dominant long range interaction []nj. Finally, the good agreement between the the- 
ory and experiment provides the first experimental confirmation of the breakdown of the 
large-quantum-number formulation of the correspondence principle |fj [| . 

I thank Anthony F. Starace and Tom Kvale for careful reading of the manuscript and 
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FIGURES 

FIG. 1. The xf functions for an attractive 1/r 6 interaction plotted vs (e^,) 1 / 3 . Solid line: 
I = 0; Dash line: I = 1; Dash-dotted line: I = 2 dotted line: I = 3. The bound spectra of any 
potential with V(r) — » —Cq/t 6 at large distances is given by the crossing points of this same set of 
functions with a set of short range parameters K c (e,l). For systems that satisfy (3§ 3> r^, K c (e,l) 
is approximately an l-independent constant in the threshold region, and the bound spectra of 
different angular momenta are given by the crossing points of xf functions with a single horizontal 
line representing K c = const. The solid horizontal line in this figure corresponds to K c = 0.3106. 
Its crossing points with Xu a l° n g with the energy scaling factor determined by Cq = 4533 a.u. and 
/i = 77392.368 a.u., give the binding energies of 85 Rb2 presented in the last columns of Table | and 
Table 0. 

FIG. 2. Partial scattering cross sections for a pair of spin-polarized 85 Rb atoms in state F\ = 3, 
Mp\ = 3, F2 = 3, and Mp2 = 3. (3§ = 162.7 a.u. is determined with n = 77392.368 a.u. and 
C 6 = 4533. 
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TABLES 

TABLE I. Comparison of energies, in GHz, of the last four bound states of 85 Rb2 characterized 
by quantum numbers F\ = 3, F2 = 3, F = 6, 1 = 2, T = 8. 



Experiment 8 



Theory I 



Theory II 



Theory IIP 




1 

2 
3 



-0.16 ±0.03 
-1.52 ±0.03 
-5.20 ±0.03 
-12.22 ±0.06 



-0.15 
-1.50 
-5.16 
-12.21 



-0.1513 
-1.520 
-5.222 
-12.37 



-0.1539 
-1.520 
-5.200 
-12.30 



a From Q. 
b From §. 

C C 6 = 4426 a.u. and K c = 0.2841. 
d C 6 = 4533 a.u. and K c = 0.3106. 

TABLE II. Bound state energies (in 

GHz) of 85 Rb 2 with quantum numbers F 1 = 3,F 2 = 3, F = 6, 1, T = F + / with / = 0,2,4,6, 
respectively. They are predicted using Cq = 4533 a.u. and K c = 0.3106. 



I = 



I = 2 



I = 4 



I = 6 



-0.2341 
-1.678 
-5.434 
-12.61 
-24.30 



-0.1539 
-1.520 
-5.200 
-12.30 
-23.91 



-1.161 
-4.659 
-11.58 
-23.01 



-0.6244 
-3.826 
-10.46 
-21.61 
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Figure 1 
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Figure 2 



